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ABSTRACT 

Aims. We used a sample of 1 1 nearby relaxed clusters of galaxies observed with the X-ray instruments XMM-Newton (EPIC) pn and 
MOS, Chandra ACIS-S and ACIS-I and BeppoSAX MECS to examine the cross-calibration of the energy dependence and normali- 
sation of the effective area of these instruments as of December 2009. We also examined the Fe XXV/XXVI line ratio temperature 
measurement method for the pn and MOS. 

Methods. We performed X-ray spectral analysis on the XMM-Newton and Chandra data for a sample of 1 1 clusters. We obtained 
the information for BeppoSAX from DeGrandi & Molendi (2002). We compared the spectroscopic results obtained with different 
instruments for the same clusters in order to examine possible systematic calibration effects between the instruments. 
Results. We did not detect any significant systematic differences between the temperatures derived in the 2-7 keV band using the 
different instruments. Also, the EPIC temperatures derived from the bremsstrahlung continuum agreed with those obtained from the 
Fe XXV/XXVI emission line ratio, implying that the energy dependence of the hard band effective area of the above instruments is 
accurately calibrated. This also indicates that deviations from ionisation equilibrium and a Maxwellian electron velocity distribution 
are negligible in the regions studied in the cluster sample. We thus consider the IACHEC sample of clusters of galaxies as standard 
candles for the calibration of the energy dependence of the hard band (2-7 keV) effective area of X-ray telescopes. On the other 
hand, the hard band EPIC/AOS fluxes disagreed by 5-10% (i.e. at 6-25cr level) which indicates a similar level of uncertainty in 
the normalisations of the effective areas of these instruments in the 2-7 keV band. In the soft energy band (0.5-2.0 keV) there are 
greater cross-calibration differences between EPIC and ACIS. We found an energy-dependent increase of ACIS versus pn bias in the 
cross-calibration of the effective area by ~ 10% in the 0.5 - 2.0 keV band. This amounts to a systematic difference of ~20% in the 
temperatures measured by the ACIS and the EPIC-pn cameras in this band. Due to the high statistical weight of the soft band data, 
the 0.5-7.0 keV band temperature measurements of clusters of galaxies with EPIC/XMM-Newton or ACIS/Chandra are uncertain by 
~ 10-15% on average. These uncertainties will also affect the analysis of the wide band continuum spectra of other types of objects 
using ACIS or EPIC. 

Key words. Instrumentation: miscellaneous - Techniques: spectroscopic - Galaxies: clusters: intracluster medium - X-rays: galaxies: 
clusters 
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f^*) ' 1 . Introduction mation on the calibration. The brightness is not high enough to 

, introduce pile-up problems which is the case with many types 

, This paper belongs to a series of studies aimed at assessing the of ca i ibration sourC es. On the other hand, the extended nature of 

■ cross-calibration status among operational and past X-ray mis- clusters of ga i axies introduces additional complexity due to e.g. 

. £h ", s i° ns promoted by the International Astronomical Consortium p 0m t spread function and vignetting 
.for High Energy Calibration (IACHEC^. In this work, we 

5^ ■ mainly explore the use of clusters of galaxies for the cross- xhe x . ray spectra of galaxy clusters C an typically be de- 

C3 ; calibration of the effective area (defined here as the product of scribe d by the combination of collisionally ionised plasma emis- 

the mirror effective area, the detector quantum efficiency and the sion componen ts, whose physics is fairly well known. While the 

filter transmission). Unlike many astronomical targets used for energy dependence (i.e. the shape) of the effective area primar- 

X-ray calibration, clusters of galaxies are stable. There is no si- jj y i m p acts the determination of the cluster temperature, the nor- 

multaneity requirement for observing clusters of galaxies with ma i isation of the effective area affects the determination of the 

different instruments. Thus, it is possible to obtain large cross- emission measure and flux. Comparing temperatures and fluxes 

calibration samples of clusters of galaxies which enables a statis- meaS ured by different instruments for the same cluster yields in- 

tically meaningful analysis of systematic effects between X-ray formation on the cross-calibration of the effective areas. For the 

instruments. Clusters of galaxies are relatively bright (L~ 10 45 hottest cmste rs, the H-like and He-like FeKa lines (Fe XXV and 

erg s ) and thus, the measurements can be performed to a high Fe X XVI) are measurable and thus the temperature-dependent 

statistical precision which translates into equally precise infer- Fe XXV/XXVI line flux ratio provides an additional diagnostic 

tool for the cross-calibration. These lines cover a narrow energy 

Send offprint requests to: J. Nevalainen band where the uncertainties of the calibration of the shape of 

1 http://web.mit.edu/iachec the effective area do not substantially affect the measurement. 
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Table 1. Information on the observations of the sample. 









XMM-Newton 










Chandra 




name 


ID 


Obs. date 


Optical filter 


exposure time (ks) 


lc filter" 


ID 


Obs. date 


exposure time (ks) 






yyyy-mm-dd 




pn 


MOS1 


MOS2 






yyyy-mm-dd 


ACIS-S 


ACIS-I 


A1795 


0097820101 


2000-06-26 


thin 


23 


37 


34 


flare 


6160 


2005-03-20 


14 


- 


A2029 


0111270201 


2002-08-25 


thin 


8 


12 


12 


flare 


4977 


2004-01-09 


78 




A2052 


0109920101 


2000-08-21 


thin 


23 


30 


30 


±20% 


5807 


2006-03-26 


127 




A2199 


0008030201 


2002-07-04 


thin 


12 


14 


14 


±20% 


497 


2000-05-13 


16 




A262 


0109980101 


2001-01-16 


thin 


16 


23 


24 


±20% 


2215 


2001-08-03 


30 




A3112 


0105660101 


2000-12-24 


medium 


17 


22 


23 


±20% 


2516 


2001-09-15 


11 




A3571 


0086950201 


2002-07-29 


medium 


7 


19 


20 


flare 


4203 


2003-07-31 


8 




A85 


0065140101 


2002-01-07 


medium 


9 


12 


12 


±20% 


904 


2000-08-19 




38 


Coma 


0153750101 


2001-12-04 


medium 


17 


21 


21 


±20% 


9714 


2008-03-20 




30 


HydraA 


0109980301 


2000-12-08 


thin 


7 


18 


19 


flare 


4969 


2004-01-13 


88 




MKW3S 


0109930101 


2000-08-22 


thin 


28 


35 


35 


±20% 


900 


2000-04-03 




56 



Notes. <n) The method for cleaning the particle flare periods: "±20%" refers to ±20% filter around the quiescence while "flare" filter denotes 
clusters where the mild flares are allowed. 



In this work, we use the clusters A1795, A2029, A2052, 
A2199, A262, A3112, A3571, A85, Coma, HydraA and 
MKW3S observed with satellites (instruments) XMM-Newton 
(pn and MOS) , Chandra (ACIS-S and ACIS-I) and BeppoSAX 
(MECS) (see Table[TJ. We compared the temperatures and fluxes 
in order to determine the status of the cross-calibration of these 
instruments for the implementation of the calibration informa- 
tion as of December 2009. 



2. The cluster properties and selected regions 

The chosen clusters are nearby (z<0.08), and bright (the flux 
in the 2-7 keV band ~ 10~ 12 - 10~ n erg s _I crrT 2 ) and most 
host a cool core (see Table 0. These clusters exhibit no major 
merger signatures in the central regions. The Galactic hydrogen 
column density is less than 6 x 10 20 cm -2 in all of the clusters. 
The temperature of the gas ranges from ~2 keV to ~10 keV in 
the cluster sample. 

The most significant temperature variation in the relaxed 
cool core clusters is due to the radial decline of the gas tem- 
perature and the central cooling. To minimise the effect of the 
radial dependence of the gas temperature between observations 
by different instruments, we aim to extract spectra within the 
same annular regions near the peak of the cluster temperature 
profile for a given cluster. This is complicated because of the 
different obscured areas due to, e.g., CCD gaps in the different 
instruments. The full annulus set by the chosen inner and outer 
extraction radii is obscured by ~15% (pn) and ~5% (MOS), on 
average, in our sample. ACIS-I observations of A85 and Coma 
cover ~85% of the full annulus, while the ACIS-I observation 
of MKW3S and all the ACIS-S observations cover 100% of the 
full annulus. We assume in this work that the differences in the 
covered areas do not produce significant effects on the derived 
gas temperatures. We minimise the emission from bright point 
sources by excluding circular regions with radii of 30" centered 
on the same sky positions in EPIC and ACIS (see Table[3]l. 

We treated the XMM-Newton EPIC instruments pn, MOS1 
and MOS2 separately, except in some relevant cases where we 
analysed the combined MOS 1 and MOS2 data for a given cluster 
and thus formed a single comparison group we identify as MOS. 
When comparing the temperatures from EPIC instruments only, 



Table 2. General information on the sample. 



name 


NH" 


z b 


R.A.(J2000) C 


Dec (J2000) 1 ' 




[10 20 cm- 2 ] 




hh mm ss 




A1795 


1.19 


0.0625 


13 48 53.0 


26 35 25.0 


A2029 


3.25 


0.0773 


15 10 56.2 


05 44 40.7 


A2052 


2.71 


0.0355 


15 16 44.5 


07 01 19.7 


A2199 


0.89 


0.0302 


16 28 38.5 


39 33 06.8 


A262 


5.67 


0.0163 


01 52 46.0 


36 09 09.1 


A3 112 


1.33 


0.0753 


03 17 57.7 


-44 14 18.3 


A3571 


4.25 


0.0391 


13 47 28.6 


-32 51 54.8 


A85 


2.78 


0.0551 


00 41 50.5 


-09 18 10.8 


Coma 


0.87 


0.0231 


12 59 35.7 


27 57 34.0 


HydraA 


4.60 


0.0539 


09 18 05.7 


-12 05 43.5 


MKW3S 


2.68 


0.0450 


15 21 51.4 


07 42 22.9 



Notes. <n) LAB weighted average (Kalberla et al., 2005). {b) From 
NASA Extragalactic Database. (r) The adopted coordinates of the clus- 
ter centre based on the location of the X-ray brightness peak in the 
XMM-Newton pn images. 



we used the data extracted from an annular region with an outer 
radius, r out , of 6 arcmitfl 

While the full sample of 1 1 clusters is observed simultane- 
ously with the XMM-Newton EPIC instruments, only a few clus- 
ters have been observed by both ACIS-S and ACIS-I. Our analy- 
sis of the ACIS-S and ACIS-I data of Coma and A1795 showed 
that the broad band temperatures are consistent within the 3% 
statistical uncertainties and that the broad band fluxes may dif- 
fer by up to 2-3%. Due to these close agreements, we combine 
the ACIS-S and ACIS-I results as a single comparison group we 
identify as ACIS. This group contains all of the 1 1 clusters. The 
small FOV of ACIS-S (~ 8' x 8') sets the limit for the outer radii 
Tout of the spectrum extraction. Depending on the placement of 
the cluster centre in the chip, the outer radii vary between 2- 
3'. Thus, for EPIC/ACIS-S comparison we used data for a given 
cluster extracted from an annular region with the outer radius set 
by ACIS-S (see Table g). 

Furthermore, for comparison with BeppoSAX MECS, we 
will use the published 2-10 keV band best-fit single temperature 
MEKAL values obtained in the central 0-2-4-6 arcmin regions 

2 The size of the inner extraction radius depends on several issues 
(e.g. PSF scatter and central cooling) and will be discussed in detail in 
Sects. [6T]and|673] 
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(deGrandi & Molendi, 2002). We refer to this paper for the de- 
tails of the data analysis, including the correction for the point 
spread function. For comparison with MECS, we extracted pn 
spectra from the same regions as used for MECS. Due to the 
limitations of the exposure time and the background, we were 
able to use 12 regions in 6 clusters in common with the pn and 
MECS (see Sect. 17721, 

3. Data processing 

For each cluster we used the data obtained by such on-axis point- 
ing (see Table Q} which was publicly available in May 2007 and 
has the largest useful exposure time after the flare-filtering (ex- 
cept that for Coma we used an ACIS-I observation from the year 
2008). 

3.1. XMM-Newton 

We processed the raw XMM-Newton data with the SASv9.0 tools 
epchain and emchain with the default parameters in order to pro- 
duce the event files. We used the latest calibration information as 
of December 2009. We also generated the simulated out-of-time 
event file, which we later used to subtract the events from the 
pn spectra registered during the readout of a CCD. We filtered 
the event files excluding bad pixels and CCD gaps. We further 
filtered the event files including only patterns 0-4 (pn) and 0- 
12 (MOS). We used the evselect-3.60.3 tool to extract spectra, 
images, and light curves, while excluding the regions contam- 
inated by bright point sources. We used the rmfgen-1.55.1 and 
arfgen- 1 .76.4 tools to produce the energy redistribution files and 
the effective area files. When running the arfgen tool, we used 
an extended source configuration and supplied an XMM-Newton 
image of the cluster in detector coordinates, binned in 0.5 arcmin 
pixels, for weighting the response. 

3.2. Chandra 

The data for each cluster were extracted from the Chandra 
archive and reprocessed with the CIAO 4.2 version of 
acis -process -.events tool. To examine the effects of recent cal- 
ibration updates, we analysed the data using CALDB 3.4, 
CALDB 4.1.1 (which included an updated HRMA effective 
area) and CALDB 4.2.0 (which includes an updated model for 
the contamination build-up on the ACIS filters). To filter out 
background flares, we generated light curves in the 2.5-7.0 keV 
and 9.0-12.0 keV energy bands using the data from the SI chip 
for ACIS-S3 observations and from the S2 chip for ACIS-I ob- 
servations. All times with background rates exceeding 20% of 
the mean were excluded from further analysis. 

4. Background spectra 

Our sample consists of nearby clusters that fill the whole FOV 
of the detectors and thus a local estimate for the background 
spectrum is not possible. 

4.1. XMM-Newton EPIC 

For the XMM-Newton EPIC analysis, we used blank sky-based 
estimates for the total sky+particle background spectra from 
Nevalainen et al. (2005). Utilising the > 10 keV (> 9.5 keV) light 
curves for the pn (MOS) we only accepted data from such peri- 
ods when the count rate is within ±20% of the quiescent level. 



However, in some cases, this led to very low exposure times. In 
order to improve the statistical precision for these clusters, we 
accepted the data accumulated during periods of mild flares, i.e. 
when the >10 keV count rate level was lower than 2.0 x the 
quiescent rate (see TableQ}. 

In order to account for the variability of the instrumental 
background due to cosmic rays, we used the sample of EPIC ex- 
posures taken with a CLOSED filter (Nevalainen et al. 2005) to 
extract particle background spectra at the same detector regions 
as used for the cluster data. We included this additional compo- 
nent in the fits after adjusting its normalisation so that the total 
background count rate prediction in the 10-14 keV (9.5-12 keV) 
band for the pn (MOS) matches that in the cluster observation^. 

We used the blank sky study of Nevalainen et al. (2005) to es- 
timate the relative uncertainty of our adopted background model 
in different energy bands and instruments (varying between 5% 
in the hard band for MOS and 20% in the soft band for the pn). 
We varied the background model by the appropriate amount and 
interpreted the change in the best-fit temperature as a systematic 
effect due to the background uncertainties. We propagated this 
uncertainty by adding it in quadrature to the statistical uncertain- 
ties of the temperature. The background uncertainty has some 
effect on the temperature measurement in the hard band (at ~ 
2% level) for the faintest clusters (A262, A2052, and MKW3S), 
where the background is the highest (10-15% of the cluster sig- 
nal). In the soft band, the background level is only a few per cent 
of that of the cluster signal, and thus the effect of the background 
uncertainty is negligible. 

In this work, we assumed that the sky background emis- 
sion is absorbed by the same amount of neutral hydrogen as 
in the cluster field. Thus, in principle, the difference between 
the hydrogen column density measured along the line-of-sight 
to the cluster and the average column density affecting the blank 
sky field pointings may bias the temperature measurements. We 
used A2052 as a conservative example of the column density 
bias, because it is the faintest object in our sample in the soft 
band (the background flux is ~6% of the cluster flux in the soft 
band). Assuming conservatively that the column density of the 
blank sky sample has the minimum value of all the fields used 
to generate the background file, i.e. 0.6 x 10 20 crrT 2 (Nevalainen 
et al., 2005), the fraction of the absorbed background flux dif- 
fers by ~15% at 0.5 keV from that obtained with the column 
density used for A2052 (2.7 x 10 20 cm -2 ). Thus, we overesti- 
mate the background model, i.e. underestimate the background- 
subtracted signal by ~1% at 0.5 keV for A2052. At higher en- 
ergies and for the brighter clusters the effect is smaller. Thus, 
the uncorrected column density variation in the blank sky and 
cluster fields has a negligible effect on our results. 

4.2. Chandra ACIS 

For ACIS, we extracted background images from the standard 
set of cti-corrected ACIS blank sky images in the Chandra 
CALDB (Markevitch et al., 2003). The exposure times in the 
background images were adjusted to produce the same 9.0- 
12.0 keV count rate as in the cluster observations. Except for 
A85, A2029 and A2199, all observations were done in very faint 
(VF) telemetry format and the VF background filtering was ap- 
plied using the CIAO tool acis-process-events. The same VF 
background screening was applied to the background data sets 
by only including events with "status=0". 



3 Due to the negligible effective area the cluster emission is insignif- 
icant at these energies 
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Table 3. Point source coordinates 



cluster 


R.A.(J2000) ; Dec (J2000) 




hh mm ss ; ° ' " 



A1795 


13 


48 


35 


1 


26 31 08.5 






















A2029 


15 


11 


06 


3 


05 41 21.7 


15 


10 


37 


3 


05 48 13.6 


15 


11 


00 


6 


05 49 22.1 


A2052 


15 


16 


44 


6 


07 05 13.8 


15 


16 


32 





06 58 51.1 


15 


16 


55 





07 01 56.2 


A2199 


16 


28 


26 


2 


39 33 52.6 


16 


28 


24 





39 33 21.1 


16 


29 


07 





39 32 39.3 


A262 


01 


52 


39 


4 


36 07 24.2 


01 


52 


39 


6 


36 10 16.5 


01 


52 


22 


3 


36 06 43.0 




01 


52 


32 





36 05 14.3 


01 


53 


07 


1 


36 09 01.2 


01 


53 


07 


5 


36 12 21.5 




01 


52 


39 


3 


36 12 32.3 


01 


52 


21 


6 


36 09 04.6 


01 


53 


02 


9 


36 12 41.0 


A3 112 


03 


18 


04 


3 


-44 13 48.3 


03 


18 


02 


4 


-44 16 43.3 


03 


17 


23 


6 


-44 15 06.7 




03 


17 


32 


3 


-44 17 34.5 


03 


17 


52 





-44 18 29.1 


03 


18 


02 


1 


-44 18 01.8 


A3571 


13 


47 


33 





-32 52 59.8 


13 


47 


29 





-32 48 34.9 


13 


47 


18 


7 


-32 49 10.1 


A85 


00 


41 


40 


5 


-09 19 47.9 


00 


41 


47 


9 


-09 20 44.0 


00 


42 


00 


3 


-09 19 30.2 




00 


41 


30 


3 


-09 15 47.1 






















Coma 


12 


59 


16 


6 


27 53 42.4 






















HydraA 


09 


17 


58 


4 


-12 04 49.6 


09 


18 


05 


4 


-12 08 00.5 


09 


18 


15 


4 


-12 05 27.0 




09 


17 


53 


5 


-12 10 37.0 


09 


17 


50 





-12 04 29.8 


09 


18 


22 


4 


-12 04 35.3 




09 


18 


23 


5 


-12 02 22.2 


09 


18 


09 


1 


-12 00 58.5 


09 


18 


19 


1 


-12 02 53.0 


MKW3S 


15 


21 


56 


2 


07 37 09.4 


15 


21 


39 


7 


07 38 33.0 


15 


21 


45 


7 


07 39 34.9 




15 


22 


12 





07 40 01.8 


15 


22 


07 


9 


07 45 22.5 


15 


21 


56 


1 


07 37 03.0 




15 


21 


52 


3 


07 36 33.4 


15 


21 


36 


3 


07 40 37.5 


15 


21 


54 


1 


07 47 03.2 



Table 4. Properties of the clusters. 



name 


I'm" 


r a 


Io,l/Io,2 fc 


r * 

i core,l 




r * 
i core,2 
it 


A* 


r 500 C 


f ■ d 

A non— isot 

% 


% 


A 1795 


1.5 


2.7 


2.7 


13.7 


0.41 


49.1 


0.67 


13 


7 


0.5 


A2029 


1.5 


2.5 


3.8 


3.0 


0.42 


33.8 


0.53 


11 


9 


0.3 


A2052 


1.7 


2.5 


8.5 


42.1 


0.97 


150.0 


0.65 


23 


5 


6 


A2199 


2.0 


2.9 


6.5 


3.3 


0.43 


54.2 


0.47 


27 


6 


6 


A262 


1.6 


2.7 


8.4 


20.0 


0.68 


100.0 


0.44 


49 


4 


5 


A3112 


1.5 


2.9 


8.3 


4.7 


0.46 


28.8 


0.50 


11 


17 


0.3 


A3571 


0.0 


2.1 


4.7 


11.3 


0.27 


122.3 


0.67 


28 


5 




A85 


1.5 


3.0 


13.3 


12.1 


0.51 


137.0 


0.7 


15 


6 


0.4 


Coma 


1.0 


5.0 




630 


0.75 






50 


8 




HydraA 


1.5 


2.7 




17.6 


0.49 






15 


9 


0.9 


MKW3S 


1.5 


2.5 




25.7 


0.44 






18 


10 


2 



Notes. (a) The adopted inner and outer radii for the annular spectrum extraction region for EPIC-AGS comparison. <b) The /3-profile parameters 
are obtained in this work, except that for Coma we adopted the values from Briel et al. (1992) (r) Estimated using the formula in Vikhlinin et al. 
(2006). {d> The contribution from outside the isothermal region (r > 0.3 rsoo) to the projected extraction region. (<l> The fraction of the emission 
within r in and r out due to the PSF scatter from the central cool region in pn. 



We estimated the effect of background modelling uncertain- 
ties (i.e. the statistical uncertainty in the 9.0-12.0 keV band used 
for normalisation and the spectral variability in the band used for 
the cluster analysis) by varying the background normalisation by 
+ 10%. We propagated in quadrature the shift in the best-fit value 
to the parameter uncertainties. 

5. Analysis methods 

We performed an X-ray spectroscopic analysis for all the data of 
our cluster sample (see Sect. [A}. The temperature measurement 
is primarily driven by the shape of the bremsstrahlung contin- 
uum for clusters with kT > 2 keV. Thus, the accuracy of the 
temperature measurement is sensitive to the accuracy of the cal- 
ibration of the energy dependence of the effective area (see Sect. 
IB. 1 .21 >. The comparison of continuum temperatures measured 
with different instruments for the same cluster yields informa- 
tion on the cross-calibration accuracy of the energy dependence 
of the effective areas. The normalisation of the effective area af- 



fects the derived emission measure and flux. Thus, the compari- 
son of cluster fluxes enables an evaluation of the uncertainties of 
the cross-calibration of the normalisation of the effective areas. 

5.1. Spectral fits 

For the spectral analysis, we modelled the emission with a 
single-temperature MEKAL model (Kaastra 1992). We used the 
metal abundances of Grevesse & Sauval (1998). We used the 
optical measurements of cluster galaxies found from the NASA 
Extragalactic Database for the cluster redshift. To model the ef- 
fect of the Galactic absorption, we applied the PHABS model. 
We fixed the column density to the value obtained from 2 1 cm 
radio observations (Kalberla et al., 2005) and used the absorption 
cross sections of Balucinska-Church & McCammon (1992). 

We binned the spectra to contain a minimum of 100 counts 
per bin. We fitted the background-subtracted (see Sect.|U) spectra 
in different energy bands: 0.5-7.0 keV (wide band) , 0.5-2.0 keV 
(soft band) and 2.0-7.0 keV (hard band). We allowed the tem- 
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perature, metal abundance and the model normalisation to vary 
in the fits, unless stated otherwise. We used the best-fit emission 
measures to determine the cluster fluxes. 

5.1 .1 . Fe XXV/XXVI line ratio method 

In the spectral analysis described above, the shape of the 
bremsstrahlung continuum dominates the temperature measure- 
ment in broad energy bands. However, the temperature depen- 
dence of the Fe ionisation fraction can also be used as an addi- 
tional, independent temperature measurement. In practise, the Fe 
XXV/XXVI line flux ratio is the most useful ionisation temper- 
ature measurement in hot clusters of galaxies. The energy reso- 
lution of EPIC and ACIS instruments is adequate for resolving 
these two lines. The lines are centred at ~6.7 and ~7.0 keV in 
the rest frame of the clusters. 

We measured the ionisation temperature by fitting the spec- 
tra in a band of 0.8 keV width centred on the mean value of 
the redshifted Fe line energies (= 6.85 keV / (1+z)) with a 
MEKAL model. To yield enough channels in this narrow band, 
we grouped the data with a requirement that each channel con- 
tains a minimum of 50 counts. 

5.2. Statistical analysis 

By comparing the spectroscopic results (temperatures (T) and 
fluxes (F) ) of our cluster sample we addressed the issue of sys- 
tematic uncertainties in the calibration of the different instru- 
ments. For each cluster, and a given pair of instruments, we com- 
puted the difference of the temperature or the flux measurements 
and its statistical uncertainty, in terms of a fraction of the aver- 
age value (fj = AT/<T> and o-f T - crAr/<T>; f^ = AF/<F> and 
<Tf F = cr AF /<F>). Using these values we calculated the weighted 
mean and its statistical uncertainty at l<x level. A systematic un- 
certainty in the calibration would tend to drive the mean of the 
distributions away from zero. We evaluate the significance of the 
deviation of the weighted mean from zero in terms of its statisti- 
cal uncertainty. 



6. Minimising the number of temperature 
components 

For optimal calibration, the spectra of the chosen targets should 
be as simple as possible. We use a single-temperature plasma 
emission model in our analysis, to make a comparison between 
different instruments. Possible multi-temperature components in 
the ICM complicate the issue, because different instruments have 
different energy dependent effective areas which will weight 
the cooler and hotter emission differently (e.g. Mazzotta et al., 
2004). Thus, even a 100% accurate calibration in all instruments 
may lead to different single-temperature measurements, if a very 
complicated temperature structure is present in the cluster. This 
is unlikely in our sample because it contains the most relaxed 
nearby clusters known in literature and we only extract spectra 
near the peak of the temperature profile. We investigate in Sect. 
16. li the effects of the central cooling and the temperature decline 
at large radii present in relaxed clusters. 

6.1. Central cooling 

We aim to minimise the effect of the central cooling in our cluster 
sample by simply excluding the emission from the cool cores in 
our analysis. Vikhlinin et al. (2005,2006) showed that in most 




Fig. 1. The solid lines show the average 3D temperature profiles 
from Vikhlinin et al. (2006) for each cool core cluster in our 
sample. The models are plotted as a function of the 3D radius 
in units of rsoo, starting from the radii corresponding to ri n when 
projected to the plane of the sky, extending to r out used for ex- 
tracting the pn/MOS spectra. The solid vertical bars indicate r out 
used for extracting the spectra for EPIC/ ACIS comparison. The 
dashed vertical lines denote the approximately isothermal region 
of 0. 1-0.3 r 500 . 

clusters the central cooling is confined within the central 70 kpc. 
At the redshifts of our cluster sample, 70 kpc corresponds to 
angular diameters of 0.8-1 .7'. We constrained the inner radius of 
our spectral extraction regions to be greater than 70 kpcQ A3571 
and Coma are exceptions to this rule, because they do not exhibit 
any central cooling. We excluded the central 1' region of Coma 
to avoid the contribution from the central galaxy NGC4874. 

6.2. Radial temperature gradient 

Most clusters feature a radial temperature profile that decreases 
at large radii, as found, e.g. from the Chandra analysis by 
Vikhlinin et al. (2005,2006). Over the 3-dimensional radial 
range of 0.1 - 0.3r<;oo the cluster temperature profiles are rather 
flat, i.e. they vary by less than 15%. We estimated the value 
of rgoo for our sample by using the T - rsoo scaling relation of 
Vikhlinin et al. (2006) and approximating the maximum tem- 
perature T pea k using our single-temperature pn fits in the 2.0- 
7.0 keV band (see Table [4). For the cool core clusters, we fur- 
ther calculated the 3D temperature profile using the model for 
the sample average in Vikhlinin et al. (2006). The calculations 
showed that, except for the nearest cluster A262, most of the ra- 
dial ranges used for the spectrum extraction in our sample are 
within the rather isothermal range of ~0. 1-0.3 rsoo (see Fig.Q}. 

In order to characterise the gas density distribution, we con- 
structed surface brightness profiles using ACIS data, except that 
for Coma we used the published results from ROSAT analy- 
sis (Briel et al., 1992). We fitted the profiles with a single-/? or 
double-/? model 



/(/?)= 7 ,iX 



1 + 



F core,! 



(-3/81+1) 



+k,2X 



1 + 



(-3ft+i) 



,(D 



4 For the additional constrain for r in , arising from the PSF, see Sect. 

E3 
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where b is the projected radius. We used the best-fit parameters 
(see Table|U) to distribute the gas density with radius. We divided 
the cluster in concentric spherical shells and assigned each shell 
with a density and temperature given by the above models. We 
then intersected the spherical shells with a hollow cylinder with 
the inner and outer radii given in Table [4] which represent the 
projected spectral extraction region, and computed the intersect- 
ing volumes of each shell. Finally, we calculated the emission 
measure of each shell at different temperatures along the line of 
sight. 

We found that 4-17% of the emission in the regions used for 
the EPIC/ AOS comparison originates from outside the nearly 
isothermal region (i.e. from 0.3-1.0 rsoo, see Table @). A3 112 
has the largest fraction (17%) of its projected emission origi- 
nating from the non-isothermal region and we thus used it as a 
conservative example on the effect of the possible temperature 
bias due to the temperature gradients in our sample. We used the 
above projected temperatures and emission measures to form a 
composite spectrum model and used it to simulate data for pn, 
MOS1 and MOS2 and ACIS-S using XSPEC. We fitted the sim- 
ulated isothermal emission (originating from 0.1-0.3 rsoo ) and 
the total emission (originating from 0.1-1.0 rgoo) with a single- 
temperature MEKAL model. The exercise showed that the ef- 
fect of including the emission from 0.3-1.0 rgoo was to reduce 
the measured temperatures by ~3%. Importantly, all instruments 
yielded the same temperature in all energy bands within ~0. 1 % 
when we included the emission from 0.3-1.0 rsoo. Thus, in our 
sample, the projected lower temperature emission due to temper- 
ature gradient does not introduce a significant bias in the derived 
temperatures between different instruments or energy bands. 

6.3. PSF scatter 

A complication when comparing Chandra and XMM-Newton 
spectra is that the point spread function (PSF) of the pn (FWHM 
~ 6") and MOS (FWHM ~ 5") (see XMM-Newton Users' 
Handboolfl) is larger than that of ACIS. The tail in the King- 
profile shape of the EPIC PSF results in 90% encircled radius 
of ~40-50", while the corresponding value for ACIS is 1-2" in 
the 0.5-7.0 keV band (see The Chandra Proposers' Observatory 
Guide©. Thus, when a bright cool core is present in a cluster, 
the tail of the PSF distribution may scatter a significant fraction 
of the central emission into the regions we examine in this work. 
This effect is greater in EPIC due to its wider PSF, and this may 
introduce a bias towards lower temperatures measured with the 
EPIC instruments. 

To minimise the effect of scattered emission from the bright 
cool core, we additionally required that the inner extraction ra- 
dius for each cluster must be greater than 1.5 arcmin (i.e. ~2 
times the 90% encircled radius) for the cool core clusters. The 
inner radius of the spectral extraction region was therefore set 
to the maximum between the cool core exclusion radius (see 
Sect. 16. U and the PSF core exclusion radius defined above ex- 
cept for the clusters without the cool core. We then estimated the 
contribution of the PSF-scattered flux within our adopted spec- 
tral extraction regions. We performed these calculations for the 
pn, which has the greatest effect. We used analytical expressions 
for the PSF profile of the pn from Ghizzardi (2002). The weak 
energy dependence of the PSF produces a negligible effect in the 
0.5-7.0 keV band and we thus report the contributions assuming 
that all photons have an energy of 2 keV. 

5 http://xmm.esa.int/extemal/xmm_user_support/documentation/uhb/ 

6 http://cxc.harvard.edu/proposer/POG/html/index.html 



We convolved the Chandra surface brightness model (see 
TableHJi from a given input radial region with the pn PSF profile 
and calculated its relative contribution to the flux within a given 
output radii for pn. The calculations showed that the contribution 
from the cool inner 70 kpc radius to the total emission projected 
within ri n and r out ranges between 0.3% and 6% (see Table|4|. 

To obtain a conservative estimate for the uncorrected PSF 
effect, we further examined A2052, which has the greatest PSF 
scattering contribution. We used the published Chandra temper- 
ature profile of A2052 (Blanton et al., 2003) to approximate the 
temperature within our projected radii r; n and r out (~3.5 keV) 
and within the central 70 kpc (~2.5 keV). Using XSPEC simu- 
lations, we then produced a composite spectrum containing two 
MEKAL components with the above temperatures in the propor- 
tion estimated by the PSF scattering calculations. We then fit- 
ted the composite spectrum with a single temperature MEKAL 
model, to approximate the effect of uncorrected PSF contribu- 
tion. The spectral analysis shows that the effect of the PSF scat- 
tering is to reduce the best-fit pn temperature by ~2%, inde- 
pendent of the band choice used in the spectral analysis. Due 
to its smaller PSF, the effect is smaller for MOS and ACIS. 
Thus, the PSF scattering may cause differences in the tempera- 
ture measurements with different instruments by a maximum of 
2%, which is negligible compared to the statistical uncertainties. 



7. Hard band temperatures 

We first restricted the spectroscopic analysis to the hard band 
(2.0-7.0 keV). This was motivated by the fact that the modelling 
of the cluster emission is less complex because the galactic ab- 
sorption does not affect the emission at these energies. Also, the 
choice of 2 ke V for the break between the soft and the hard bands 
allows us to have some control on the calibration of the different 
instrumental components. The EPIC optical filter transmission 
is essentially 100% in the hard band, while it is significantly 
smaller in the soft band. The molecular contamination on the 
ACIS filters (see below) also affects only the lowest energies. 
The quantum efficiency of the pn is essentially 100% in the hard 
band, and thus the effective area of the pn primarily depends 
on the effective area of the mirrors. For MOS and ACIS, the 
quantum efficiency is less than 100% in the hard band and the 
different components cannot be distinguished. 

7.1. EPIC temperatures 

The single temperature MEKAL hard band fits to the pn, MOS 1 
and MOS2 data revealed that the pn and MOS 1 values differ by 
only ~1% on average (see Fig. [2] and Table [5]). There is an in- 
dication that the MOS2 temperatures are systematically lower 
(by ~4% and ~5%) than those obtained with the pn and MOS1, 
respectively. However, these differences are significant only at a 
modest level of 2.6<x and 3.4cr. The combined MOS data (MOS 1 
+ MOS2) yields ~2% lower hard band temperatures than the pn 
on average and the values agree within the uncertainties at the 
1 cr confidence level (see Fig. |3j. Thus, the cross-calibration of 
the energy dependence of the effective area between the pn and 
MOS is consistent within the statistical precision of these mea- 
surements, which corresponds to less than ~5% discrepancy in 
temperature. We will use the pn data in the following discussion 
for comparisons with the hard band temperatures obtained with 
other satellites, due to the smaller statistical uncertainties with 
respect to the MOS cameras. 
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2468 10 2468 10 2468 10 

pn hard band T [keV] pn hard band T [keV] M0S1 hard band T [keV] 



Fig. 2. The upper panels show the best-fit temperatures (asterisks) and the Icr uncertainties (boxes) in the hard band obtained with 
the XMM-Newton instruments pn (left panel) , MOS1 (middle panel) and MOS2 (right panel). The lower panels show the relative 
temperature difference, fj = AT/<T> (diamonds), and its l<x uncertainty for MOSl-pn (left panel), MOS2-pn (middle panel) and 
MOS2-MOS1 (right panel). The dotted and dashed lines show the weighted mean of fj + the error of the mean. 




pn hard band T [keV] 



Fig. 3. The best-fit temperatures (asterisks) and the 1 <x uncertainties (boxes) in the hard band obtained with the pn v.s. those obtained 
with MOS (i.e. combined MOS1 and MOS2 spectra). 



7.2. XMM-Newton v.s. BeppoSAX temperatures 

We found that the pn hard band temperatures only differ by ~ 1 % 
(0.9cr) from the BeppoSAX MECS values on average (see Fig. 
[4] and Table [6}. The pn values are not systematically higher or 
lower than the MECS values. 



7.3. XMM-Newton v.s. Chandra temperatures 

Our spect ral fi ts to the pn and ACIS data using CALDB 4.2.0 
(see Figs. lA.ll and IA.3b showed that the ACIS hard band tem- 
peratures differ from the pn values only by ~l% (0.6cr) on aver- 
age (see Fig.|5]and Table IA. 11 implying similar accuracy for the 
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Table 5. Best-fit hard band temperatures and lcr uncertainties 
for EPIC instruments. 



Table 6. Best-fit hard band temperatures and lcr uncertainties 
forpnandMECS. 









MOS1 


MOS2 


MOS 


pn 


name 


r in 


r ou , 


Thard 


Than] 


Thard 


Thard 




[•] 


['] 


[keV] 


[keV] 


[keV] 


[keV] 


A 1795 


1.5 


6.0 


7.0(6.8-7.2] 


6.3(6.1-6.5] 


6.6(6.5-6.8] 


6.7(6.5-6.9] 


A2029 


1.5 


6.0 


9.2(8.7-9.7] 


8.4(7.9-8.9] 


8.6(8.2-9.0] 


8.0(7.6-8.4] 


A2052 


1.7 


6.0 


3.4[3. 3-3.5] 


3.4(3.3-3.5] 


3.4(3.3-3.5] 


3.5(3.4-3.6] 


A2199 


2.0 


6.0 


4.7[4.5-4.8] 


4.5(4.4-4.6] 


4.6(4.5-4.7] 


4.7(4.6^1.8] 


A262 


1.6 


6.0 


2.5(2.4-2.6] 


2.5(2.4-2.5] 


2.5(2.4-2.5] 


2.5(2.4-2.6] 


A3112 


1.5 


6.0 


6.1|5.7-6.5| 


5.2(4.9-5.5] 


5.6(5.3-5.9] 


5.7(5.4-6.1] 


A3571 


0.0 


6.0 


7.4(7.2-7.6] 


7.0(6.8-7.2] 


7.2(7.1-7.4] 


7.4(7.2-7.6] 


A85 


1.5 


6.0 


6.5(6.2-6.8] 


7.1(6.7-7.4] 


6.7(6.4-7.0] 


7.0(6.7-7.3] 


Coma 


1.0 


6.0 


9.4(9.1-9.8] 


8.9(8.6-9.2] 


9.2(8.9-9.4] 


9.0(8.8-9.3] 


HydiaA 


1.5 


6.0 


4.8(4.5-5.0] 


4.3(4.1^1.5] 


4.5(4.3-4.7] 


5.0(4.7-5.4] 


MKW3S 


1.5 


6.0 


4.0(3.9-4.1] 


3.7(3.6-3.8] 


3.8(3.7-3.9] 


4.0(3.8-4.1] 



Notes. Spectra were extracted from concentric annuli of inner and outer 
radii given as r in and r out . The data were fitted with a single-temperature 
MEKAL model in the hard band (2-7 keV) with the absorption fixed at 
the Galactic value and the abundance allowed to vary. 




pn hard band T (keV) 

Fig. 4. The upper panel shows the best-fit temperatures (aster- 
isks) and lcr uncertainties (boxes) in the hard band obtained 
with the XMM-Newton pn and BeppoSAX MECS instruments. 
The lower panel shows the relative MECS-pn temperature dif- 
ference, fj = AT/<T> (diamonds), and its 1 cr uncertainty. The 
dotted and dashed lines show the weighted mean of fj + the 
error of the mean. 



cross-calibration of the energy dependence of the effective area 
between ACIS and EPIC in the hard band. 



8. Fe XXV/XXVI line ratio 

In this Chapter, we explore the Fe XXV/XXVI line ratio method 
for the determination of gas temperatures. 



8.1. Statistical issues 

We obtain the Fe XXV/XXVI line ratio temperatures by fit- 
ting the data with the MEKAL model in the ~6-7 keV energy 









pn 


MECS 


name 


I'm 


Toil l 


Thard 


Thard 




['] 


['] 


[keV] 


[keV] 


A1795 


2.0 


4.0 


6.6[6.4-6.9] 


6.2[6. 1-6.4] 


A1795 


4.0 


6.0 


7.0[6.4-7.6] 


6.3[6.0-6.5] 


A2029 


2.0 


4.0 


8.0[7.5-8.5] 


8.0[7.7-8.4] 


A2199 


2.0 


4.0 


4.8[4.7-5.0] 


4.5[4.5-4.6] 


A2199 


4.0 


6.0 


4.7[4.5-4.9] 


4.6[4.5-4.7] 


A3571 


0.0 


2.0 


7.5[7. 1-7.8] 


7.5[7.3-7.8] 


A3571 


2.0 


4.0 


7.7[7.3-8.1] 


7.6[7.4-7.8] 


A3571 


4.0 


6.0 


7.0[6.6-7.5] 


7.7[7.4-8.0] 


A85 


2.0 


4.0 


7.5[7.0-8.1] 


6.9[6.7-7.1] 


Coma 


0.0 


2.0 


9.7[9.1-10.4] 


9.2[8.9-9.7] 


Coma 


2.0 


4.0 


9.3[9.0-9.7] 


10.0[9.6-10.4] 


Coma 


4.0 


6.0 


8.9[8.6-9.2] 


9.2[9.0-9.5] 



Notes. The values are obtained using data extracted from concentric 
annuli of inner and outer radii given as r; n and r m „. The data were fitted 
with a single-temperature MEKAL model in the hard band (pn: 2-7 
keV, MECS:2-10 keV). 

Table 7. Systematic effects on the temperature measurement. 



source 


cr T /T 


comments 




MEKAL" 


background 


<2% 


EPIC and ACIS, propagated 


PSF 


« 2% 


EPIC and ACIS, not propagated 


T gradient 


<0.1% 


EPIC and ACIS, not propagated 


S 


< 2% 






Fe XXV/XXVI fc 


limited counts 


<2% 


EPIC 


energy resolution 


2% 


pn only 


gain 


1-3% 


EPIC 


MEKAL/APEC 


2% 


EPIC 


S 


4% 


propagated 



Notes. <D) Temperature measurements using the MEKAL model in the 
hard band (6) Temperature measurements using the Fe XXV/XXVI line 
ratio 



band, depending on the cluster redshift (as explained in Sect. 
15.1.1b . We first attempted to estimate whether the accuracy of 
the method depends on the number of photons available for the 
analysis. For this, we performed a spectral analysis of data simu- 
lated with the MEKAL model in XSPEC. We used the pn energy 
redistribution and auxiliary files of A1795 for the simulations. 
We used a range of 5-10 keV for the temperatures in the simula- 
tions and kept the metal abundance at 0.3 Solar and a redshift of 
zero in the simulations. We varied the number of input counts in 
a 6.45-7.25 keV energy band by scaling the input model accord- 
ingly. For each temperature and number of counts, we performed 
1000 simulations, including the scatter due to counting statistics. 

We binned the simulated spectra to contain a minimum of 
50 counts in each channel, as we do when fitting the actual data. 
We then fitted each of the simulated spectra in the 6.45-7.25 keV 
band and thus obtained a distribution of best-fit temperatures for 
each input temperature and number of counts. We used these dis- 
tributions to determine the median temperatures and compared 
them with the input values. The comparison showed that while 
the derived temperatures agreed with the input values within 1 % 
at the highest numbers of counts (> 3000), at the low count num- 



8 



J. Nevalainen et al.: IACHEC cluster calibration 




pn hord bond T [keV] 

Fig. 5. The best-fit temperatures (asterisks) and the lcr uncertainties (boxes) in the hard band obtained with the XMM-Newton pn 
detector v.s. those obtained with the Chandra ACIS detector using CALDB 4.2.0. 



bers there is a tendency towards higher temperatures, especially 
for the lowest input temperatures (see Fig. |6]). This effect was 
identical regardless of whether we used^ 2 or C-statistics. At the 
lowest temperatures, the Fe XXVI line emission is quite weak 
and it is possible that the Fe XXVI line-like features are created 
at random due to low statistics. These artificial features are then 
interpreted as real Fe XXVI emission and the temperatures are 
biased high, as measured above. 

We adopt as a selection criterium for our sample that the 
uncertainty due to the limited number of counts biases the tem- 
perature measurements by less than 2%. We consider this 2% 
uncertainty below when discussing the various sources of uncer- 
tainty in the temperature measurements. For clusters with tem- 
perature above ~7 keV, this requirement imposes a lower limit 
of ~1000 counts in the 6.45-7.25 keV band. For cooler clusters, 
the requirement for the minimum number of counts is higher. 

These limits are quite severe when compared to the observed 
numbers of counts in our sample. This is especially true for the 
EPIC-AGS comparison where the inner extraction radius is lim- 
ited due to our minimisation of PSF scattering, and the outer 
extraction radius is limited by the ACIS FOV. Consequently, 
the number of counts in the ~6-7 keV band in most clusters 
is well below the above limits. Thus, with the current exposures 
we cannot perform a statistically meaningful comparison of Fe 
XXV/XXVI temperatures between EPIC and ACIS. On the other 
hand, due to the larger regions used in the pn-MOS comparison, 
the number of counts in the ~6-7 keV band exceeds 1000 in 
most clusters. Out of these clusters, we chose those with kT> 5 
keV for our analysis, i.e. A1795, A2029, A3112, A3571, A85 
and Coma. 
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x 1000 counts in the 6.45-7.25 keV band 



Fig. 6. The relative bias of the temperature derived by fitting the 
simulated spectra in the 6.45-7.25 keV band as a function of 
the number of counts for the input values kT=5-7 keV (dotted 
blue lines) and for kT=8-10 keV (solid lines). The dash-dot line 
shows the allowed upper limit of 2% for the bias in the sample. 

8.2. Systematic uncertainties 

The accuracy of the Fe XXV/XXVI line ratio measurement ad- 
ditionally depends on the accuracy of the absolute energy scale 
and the energy resolution used to generate the redistribution ma- 
trix in the ~6-7 keV band. Also, uncertainties in the redshift 
measurement and the atom physics used to generate the MEKAL 
model have some effect to the interpreted temperature for a given 
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Fig. 7. The best-fit MEKAL models (solid line) in the ~6-7 ke V band data which contain the Fe XXV and Fe XXVI emission lines 
and statistical uncertainties (crosses). The pn (left panel) and MOS (right panel) data and the models are arbitrarily scaled for the 
clarity of the display. 



line ratio data. We found that these effects (see Appendix [B] 
for the details), combined with the bias due to limited number 
of counts (Sect. I8.U may affect the EPIC temperatures at the 
~4% level (see Table |7J. We will consider this uncertainty and 
propagate it along with the statistical uncertainties of EPIC Fe 
XXV/XXVI temperatures when appropriate in the following. 

8.3. Results 

We then fitted the EPIC data of the subsample described in Sect. 
I8.1l (see Fig. [7]) in the ~6-7 keV band with a single-temperature 
MEKAL model. To obtain a similar level of counts as in the pn 
instrument, we co-added the MOS 1 and MOS2 data for compar- 
ison. 

The relatively small number of photons in the narrow energy 
band (6.45-7.25 keV in the rest frame of a cluster) inevitably 
limits the statistical accuracy of the derived temperatures. The 
relative lcr statistical uncertainty on the temperatures is at most 
20% (see Table©. In the 6.45-7.25 keV energy band, the metal 
abundance and the emission measure are highly degenerate: a 
lower normalisation for the continuum is approximately com- 
pensated with a higher metal abundance, since the continuum 
cannot be well separated from the line profile in this narrow en- 
ergy band (see Fig. [8]). 

On the other hand, we have already constrained the emission 
measure when fitting the hard band (2-7 keV) (Sect. [7]), typically 
within 1% at lcr level. We thus experimented by constraining 
the emission measure to that found in the hard band fit, when 
fitting the ~6-7 keV band. The resulting statistical uncertainties 
for the temperature are 3-10%, instead of 5-20% in the case of a 
free emission measure (see Table|8]l. Thus, with this method, we 



reach the level of the systematic uncertainties in some cases (see 
Table |7). In both cases, the temperatures are consistent within 
the statistical lcr uncertainties. Importantly, the temperatures ob- 
tained with the free or constrained emission measure in the ~6- 
7 keV band do not indicate a systematic difference (see Table 
[8}. This proves that the use of the hard band emission measure 
constrain does not introduce any bias to the EPIC temperature 
measurements. 

The pn and MOS temperatures obtained with the hard band 
emission measure constrain mostly agree within the statistical 
lcr uncertainties (see Tableland Fig . [9j ■ The values differ some- 
what for A1795 and A3571, which have the smallest statistical 
uncertainties. Considering the systematic uncertainties, the val- 
ues agree. There is no systematic difference between the derived 
temperatures for the full sample obtained with the two instru- 
ments. Thus, in the following we will refer to the Fe XXV/XXVI 
temperatures obtained with the emission measure constrain. 

8.4. Fe XXV/XXVI line ratio v.s. continuum 

The ~6-7 keV band temperatures above are given by the tem- 
perature dependence of the ionisation fraction, i.e. by the Fe 
XXV/XXVI line ratio. On the other hand, the hard band tem- 
peratures are primarily determined by the temperature depen- 
dence of the bremsstrahlung continuum. There are several rea- 
sons why the two measurements may not agree. Due to cluster 
mergers, the intracluster material may be in a non-equilibrium 
ionisation state which would change the Fe XXV/XXVI line 
ratio from that given by ionisation equilibrium assumed in 
our work (e.g. Prokhorov, 2010; Akahori & Yoshikawa, 2010). 
The mergers may also accelerate a fraction of electrons into 
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Table 8. Best-fit temperatures and statistical uncertainties for the EPIC instruments in the ~6-7 keV band. 













pn 






MOS 










EM free" 




EM constr. 


b 


EM free" 




EM constr. 


b 


name 


I'in 


Tout 


T 


cr r /T 


T 


a T [Y 


T 


crr/T 


T 


cr T fT 




['] 


['] 


[keV] 




[keV] 




[keV] 




[keV] 




A1795 


1.5 


6.0 


6.1 [5.7-6.5] 


7% 


6.6[6.4-6.8] 


3% 


5.7[5.3-6.1] 


7% 


6.1 [5.9-6.3] 


3% 


A2029 


1.5 


6.0 


7.7[7. 1-8.4] 


9% 


7.7[7.3-8.2] 


6% 


9.0[8.3-9.7] 


8% 


8.4[7.8-8.9] 


7% 


A3112 


1.5 


6.0 


5.6[4.3-6.7] 


21% 


5.6[5.2-6.2] 


9% 


7.2[6.4-8.1] 


11% 


6.4[5.9-6.9] 


&% 


A3571 


0.0 


6.0 


7.5[7.0-7.9] 


6% 


7.5[7.2-7.9] 


4% 


6.8[6.5-7.1] 


5% 


6.8[6.6-7.0] 


3% 


A85 


1.5 


6.0 


7.4[6.6-8.1] 


10% 


7.3[6.8-7.7] 


7% 


6.4[5.4-7.2] 


15% 


6.4[5.8-6.8] 


8% 


Coma 


1.0 


6.0 


9.3[8.9-9.8] 


5% 


9.0[8.7-9.3] 


4% 


9.2[8.6-9.7] 


6% 


8.9[8.5-9.3] 


4% 



Notes. The uncertainties are given at lcr. The temperatures are derived using data extracted from concentric annuli of inner and outer radii given 
as r in and r out . The data were fitted with a single-temperature MEKAL model in the ~6-7 keV band allowing the emission measure to be free or 
ib) constraining the emission measure within that derived in the hard band. 



suprathermal velocities, which will produce deviations from a 
Maxwellian electron velocity distribution assumed when mod- 
elling the continuum emission with a bremsstrahlung model 
(Prokhorov et al., 2009). Also, the non-Maxwellian electron dis- 
tribution will affect the Fe XXV/XXVI line ratio (e.g. Kaastra 
et al., 2009). Furthermore, a relativistic electron population cre- 
ated by a strong merger shock may produce an additional contin- 
uum emission component via inverse compton scatter of cosmic 
microwave background photons (e.g. Sarazin, 1988). These ef- 
fects are probably small in our sample because we used the most 
dynamically relaxed clusters. Also, we excluded the cluster out- 
skirts, which may still be accreting material and where ionisation 
equilibrium may not yet have been established. 

With the current XMM-Newton EPIC data we are able to test 
the above possibilities by a direct comparison of the tempera- 
tures derived from both the continuum and the Fe XXV/XXVI 
line ratio. To keep the two temperature measurements as inde- 
pendent as possible, we derived the continuum temperature us- 
ing a MEKAL fit in the 2.0-6.0 keV band. In this band, the line 
emission is negligible for hot clusters and thus the metal abun- 
dance may be inaccurately determined. We examined this issue 
by fixing the metal abundance to values in the range 0.0-0.5 so- 
lar. Examining the best-fit models showed that the temperatures 
varied by less than 1 % for a given cluster when varying the abun- 
dance. Thus, uncertainties in the metal abundance do not signif- 
icantly bias the temperature measurements. 

The comparison of the temperatures derived by fitting the 
2.0-6.0 keV and ~6-7 keV band reveal a good agreement be- 
tween the two methods (see Fig. ITOb : the PN temperatures dif- 
fer by -1% (0.2cr) while MOS differences are -7% (1.6cr) on 
average. There is no systematic difference between the measure- 
ments. This implies that the calibration of the energy dependence 
of the effective area of the EPIC instruments in the hard band is 
accurate and that the deviations from the ionisation equilibrium 
state and Maxwellian electron velocity distribution in these clus- 
ters are negligible. Given our evidence for the Fe XXV/XXVI 
line ratio measurement being insensitive to calibration uncertain- 
ties in the effective area (see Appendix IB. U . the above finding 
could be useful for the calibration of the hard band of future X- 
ray missions: the hard band calibration should be fine-tuned to 
yield the same temperature as does the Fe XXV/XXVI line ratio. 

9. Soft band temperatures 

We then performed a spectroscopic analysis of the data in the 
soft band (0.5-2.0 keV). We used a single-temperature MEKAL 
model, as in the analysis of the hard band above. We allowed 




n 1 1 1 6.5 r ■ 

2.6 2.8 3.0 3.2 3.4 3.6 3.8 2.6 2.8 3.0 3.2 3.4 3.6 3.8 

MEKAL normalisation x 100 MEKAL normalisation x 100 



Fig. 8. The 2-parameter confidence contours at lcr for A2029 
obtained with the MEKAL fits to the pn data in the -6-7 keV 
band. 




5 6 7 8 9 

pn Fe XXV/XXVI rotio T [keV] 



Fig. 9. The best-fit temperatures (asterisks) obtained with the 
MEKAL fits in the ~6-7 keV band data of the pn v.s. those ob- 
tained with MOS. The solid boxes indicate the statistical uncer- 
tainties at lcr level while the dotted boxes include the systematic 
uncertainties. 



the temperature, metal abundance and emission measure to vary 
independently from the hard band values presented above. This 
was required in order to obtain statistically acceptable fits in the 
soft band (see Figs. I A. 2 1 and \AA\ for pn and ACIS). Note that the 
formally accurate models for the soft and hard bands, which are 
adequate for calibration work, do not yield a physically consis- 
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5 6 7 8 9 

Fe XXV/XXVI line ratio T (keV) 



5 6 7 8 9 

Fe XXV/XXVI line ratio T (keV) 



Fig. 10. The upper panels show the best-fit temperatures (aster- 
isks) and lcr uncertainties (boxes) for the Fe XXV/XXVI line 
ratio temperatures v.s. the 2-6 keV band continuum for the pn 
(left panel) and MOS (right panel). The lower panels show the 
relative temperature difference fj (i.e. the 2-6 keV band contin- 
uum temperature minus the Fe XXV/XXVI line ratio tempera- 
ture, divided by the mean of the two) and its lcr uncertainty for 
the pn (left panel) and MOS (right panel). The dotted and dashed 
lines show the weighted mean of fj + the error of the mean. 



tent modelling of the full energy band. In order to keep the mod- 
elling and the temperature comparison simple, we excluded the 
coolest cluster A262 from this comparison, because its strong 
line emission at ~1 keV would require multi-component mod- 
elling for an accurate description of the data. 



9.1. pn v.s. MOS 

We found that while the pn and MOS2 agree very well (average 
difference is -1% , i.e. 0.5(f), the MOS1 values are -5% (-3.7 
cr) higher than those of the pn and MOS1 (see Fig. ITTb . These 
differences indicate some remaining cross-calibration uncertain- 
ties in the EPIC instruments in the soft band which we will not 
examine in more detail here. 



9.2. pn v.s. ACIS 

We found that the ACIS soft band temperatures are ~ 18% higher 
than those obtained with the pn (see Fig. [T2l . This difference 
is statistically very significant, 8.6<x, and thus there are signif- 
icant remaining ACIS/pn cross-calibration uncertainties in the 
soft band. 

To examine this difference in more detail, we convolved the 
best-fit soft band pn models through the ACIS responses and 
compared this prediction with the ACIS data. In the case of A85 
and Coma the data are spread across several CCD chips and for 
these clusters we co-added the data and formed a single average 
file for the responses and background. We scaled the models to 
match the ACIS data at 2 keV, since we are primarily interested 
in the energy dependence of the effective area here, whose cross- 
calibration is very accurate in the 2-7 keV band (see Sect. [7J. 

We found that the ratio of the ACIS data to the pn model 
prediction exhibits a systematic trend in most clusters. Fitting the 
data-to-model ratio with a linear function we found that while 




pn soft band T [keV] 




pn soft band T [keV] 



MOST soft band T [keV] 



Fig. 11. The upper panels show the best-fit temperatures (aster- 
isks) and statistical lcr uncertainties (boxes) in the soft band with 
the XMM-Newton instruments pn/MOSl (left panel), pn/MOS2 
(middle panel) and MOS1/MOS2 (right panel). The lower pan- 
els show the relative temperature difference fj ( = AT/<T>, dia- 
monds) and its 1 cr uncertainty for MOS 1 -pn (left panel), MOS2- 
pn (middle panel) and MOS2-MOS1 (right panel). The dotted 
and dashed lines show the weighted mean of fx ± the error of 
the mean. 
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pn soft band T (keV) ACIS hard band T (keV) 

Fig. 12. The upper panels show the best-fit temperatures (aster- 
isks) and lcr uncertainties (boxes) in the soft band for the pn and 
ACIS (left panel) and in the hard and soft bands for ACIS (right 
panel). The lower panels show the relative temperature differ- 
ence, fx = AT/<T> (diamonds) and its statistical lcr uncertainty 
for ACIS soft band - pn soft band (left panel) and ACIS soft 
band - ACIS hard band (right panel). The dotted and dashed lines 
show the weighted mean of f T + the error of the mean. 



the ratio is 1.0 at 2 keV by definition, it decreases to ~0.9 at 
0.5 keV in most clusters (see Fig. [T3l . This implies that either 
the effective area of the pn is underestimated, or that of ACIS is 
overestimated at 0.5 keV by ~10%. 



10. Wide band temperatures 

Since both the intrinsic cluster flux and the effective area of EPIC 
and ACIS are higher in the soft band, the statistical weight of the 
soft band data is much greater than that of the hard band. As 
shown above, the greatest differences between the accuracy of 
the EPIC and ACIS effective areas are below 2 keV which thus 
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Fig. 13. The ratio of the ACIS data to the pn model prediction using the best-fit pn model and the ACIS responses. The model has 
been normalised so that the ratio equals unity at 2.0 keV. The dashed line shows the best-fit linear relation. 



Table 9. The soft band temperatures and lcr uncertainties using 
EPIC instruments. 









pn 


MOS1 


MOS2 


name 


I'in 


1'otil 


Tsoft 


T so ft 


Tsoft 




['] 


['] 


[keV] 


[keV] 


[keV] 


A1795 


1.5 


6.0 


4.4[4.3^1.5] 


4.6[4.4^t.7] 


4.1 [4.0-4.2] 


A2029 


1.5 


6.0 


5.2[5.0-5.4] 


5.7[5.4-6.1] 


4.9[4.6-5.1] 


A2052 


1.7 


6.0 


2.5 [2.4-2.5] 


2.7[2.6-2.8] 


2.6[2.5-2.6] 


A2199 


2.0 


6.0 


3.6[3.5-3.7] 


4.0[3.8-4.1] 


3.5[3.4-3.7] 


A3112 


1.5 


6.0 


3.9[3.7^U] 


3.8[3.6-4.0] 


3.6[3.5-3.8] 


A3571 


0.0 


6.0 


5.4[5.3-5.6] 


5.5[5.3-5.7] 


5.3[5.2-5.5] 


A85 


1.5 


6.0 


4.9[4.7-5.1] 


4.9[4.7-5.2] 


4.6[4.4-4.9] 


Coma 


1.0 


6.0 


6.1 [5.9-6.2] 


6.3[6.0-6.5] 


6.3[6. 1-6.6] 


HydraA 


1.5 


6.0 


3.0[2.9-3.1] 


3.0[2.9-3.1] 


2.9[2.8-3.0] 


MKW3S 


1.5 


6.0 


2.7[2.7-2.8] 


2.9[2.9-3.0] 


2.8[2.8-2.9] 



Notes. The data are extracted from concentric annuli of inner and outer 
radii given as rj„ and r ou t- The data were fitted with a single- temperature 
MEKAL model in the soft band (0.5-2.0 keV). 



affect the temperatures derived in the wide band (0.5-7.0 keV). 
For the interest of the general user, we also evaluated the effect of 
the cross-calibration uncertainties when fitting the 0.5-7.0 keV 
band with a single-temperature model as is often done in the 
literature. Since the hard and soft band temperatures are incon- 
sistent in all clusters and for all instruments, the application of 



a single-temperature model in the wide band will result in some 
residuals. These residuals are usually hidden by the larger statis- 
tical uncertainties in detailed temperature profile measurements 
presented in the literature. 

As expected, due to the above findings (Sect. 19.21 1. ACIS 
yields systematically higher wide band temperatures than EPIC. 
On average, the ACIS/MOS1, ACIS/MOS2 and ACIS/pn differ- 
ences are ~8%, ~15% and ~14% which corresponds to a very 
high detection significance of 10. Oct 18.9<t and 18. Oct. The max- 
imum difference is ~25% (see Fig. fT4l . These values may be 
taken as a systematic uncertainty in the wide band temperature 
measurement when using EPIC or ACIS data. These uncertain- 
ties also apply to the analysis of EPIC and ACIS wide band con- 
tinuum spectra of other types of objects. Evaluation of this effect 
remains to be carried out by analysis of e.g. blazars observed 
simultaneously with XMM-Newton and Chandra . 

11. Hard band v.s. soft band temperatures 

With our data set we can also compare the temperature measure- 
ments in the hard and soft bands. In principle, this could be an 
useful tool for discovering emission components in addition to 
the bulk of the isothermal emission. 

We found that all of the XMM-Newton soft band tempera- 
tures are significantly smaller than the hard band values (see Fig. 
[I5]and Tables|5]and|2]i: The pn, MOS 1 and MOS2 soft band tem- 
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2 4 6 8 10 
pn hard band T [keV] 



2 4 6 8 10 
M0S1 hard band T [keV] 



2 4 6 
M0S2 hard band 



8 10 
T [keV] 



Fig. 15. The upper panels show the best-fit temperatures (asterisks) and statistical lcr uncertainties (boxes) in the soft and hard bands 
for XMM-Newton instruments pn (left panel), MOS1 (middle panel) and MOS2 (right panel). The lower panels show the relative 
soft-hard temperature difference fr( = AT/<T>, diamonds) and its statistical lcr uncertainty for pn (left panel), MOS1 (middle 
panel) and MOS2 (right panel). The dotted and dashed lines show the weighted mean of f T + the error of the mean. 



l — i — i — i — i — i — i — i — i — i — i — i — r- 




3.0 4.0 5.0 6.0 7.0 8.0 9.0 
pn wide bond T [keV] 

Fig. 14. The upper panels show the best-fit temperatures (aster- 
isks) and lcr uncertainties (boxes) in the wide band using pn or 
ACIS data. The lower panel shows the relative ACIS-pn tem- 
perature difference fj ( = AT/<T>, diamonds) and its statistical 
lcr uncertainty. The dotted and dashed lines show the weighted 
mean of fj + the error of the mean. 



peratures are ~30-35% (~22-30cr) lower than those obtained in 
the hard band. Also, most of the ACIS soft band temperatures are 
lower than the hard band values, but by a smaller amount, 12% 
(1.6cr, see Fig.[T2land Table lA.ll . It would be useful to examine 
whether this effect is present in the data of objects of different 
type. While the scatter is quite large and no clear patterns can 
be seen, there is some indication that the relative difference in- 
creases with temperature. 



These systematic differences between the hard and soft band 
temperatures in all instruments could be explained by devia- 
tions from the isothermal emission. Since we have minimised 
the gas cooling effect, among the candidates for the excess 
emission are multi-temperature gas due to, e.g., clumping or 
cluster mergers and the Inverse Compton scatter of Cosmic 
Microwave Background photons from relativistic cluster elec- 
trons (see Durret et al. 2008 and Rephaeli et al. 2008 for recent 
reviews). Typically this excess emission flux is ~10% of the 
isothermal component. We found above that there are remain- 
ing uncertainties in the calibration of the effective area in the 
studied instruments at this level in the soft band. Due to these 
possible astrophysical effects in the clusters of galaxies, they are 
not a good choice for calibrating the 0.5-2.0 keV band and other 
sources should be used (e.g. white dwarfs and isolated neutron 
stars). 



12. Fluxes 

We then examined the fluxes given by the best-fit MEKAL mod- 
els in the hard and soft bands. We calculated the statistical lcr 
uncertainties of the flux measurements by considering the vari- 
ation of all the free MEKAL parameters. We limited the pn- 
MOS comparison within the smaller outer extraction radii used 
for ACIS (see Sect. 0. Due to the CCD gaps and bad columns 
in the central regions the fraction of the useful detector area to 
the full annulus given by the extraction radii is less than 100% 
(except for ACIS-S), and varies between the instruments for a 
given cluster. Also the exclusion of the point source regions re- 
duces the useful detector area. To enable a meaningful compar- 
ison of the fluxes derived with different instruments, we divided 
the measured flux by the fraction of the covered full annulus to 
recover the intrinsic cluster flux measured by each instrument. 
Since the emission is not constant with radius, the linear scaling 
may produce some differences for the fluxes derived with differ- 
ent instruments for a given cluster due to the differences in the 
covered areas (see below). 
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Fig. 16. The relative flux difference, fp - AF/<F> (diamonds) and its statistical uncertainties in the hard band (upper row) and in 
the soft band (lower row) are shown as a function of pn hard band temperature for different instrument pairs. The dotted and dashed 
lines show the weighted mean of fp + the error of the mean. 



Since the half-width of the extraction region is in most cases 
comparable to the 90% encircled radius (~40-50 arcsec) of 
EPIC we expect that ~10% of the flux originating within the 
extraction region will end up outside the studied region (= out- 
flux). We confirmed this by convolving the surface brightness 
models with the pn PSF and calculating the number of counts 
in different regions. However, the flux originating outside of the 
extraction region and ending up to the extraction region (= in- 
flux) compensates very accurately the out-flux. Thus the net ef- 
fect is that the true cluster flux within our extraction regions is 
affected by less than 0.5% due to PSF scatter in our sample. The 
very accurace cancellation of the in-flux and out-flux is a con- 
sequence of an approximately constant count rate within annuli 
of constant width towards smaller or bigger radii within the PSF 
radius of 1 arcmin. This in turn is due to the opposite effects of 
decreasing surface brightness and increasing area of the annuli 
with an increasing radius. 



Table 10. The fluxes in the hard band (2-7 keV) given by the 
best-fit single-temperature MEKAL model. 









pn 


MOS1 


MOS2 


ACIS 


name 


•'in 
['] 


r ou , 
['] 


flux" 


flux" 


flux" 


flux" 


A 1795 


1.5 


2.7 


1.08±0.01 


1.08+0.01 


1.07±0.01 


1.16±0.02 


A2029 


1.5 


2.5 


1.14±0.02 


1.21±0.02 


1.14±0.02 


1.27±0.01 


A2052 


1.7 


2.5 


0.31±0.01 


0.33±0.01 


0.31±0.01 


0.34±0.01 


A2199 


2.0 


2.9 


0.86±0.01 


0.90±0.01 


0.89±0.01 


0.89±0.01 


A262 


1.6 


2.7 


0.21±0.01 


0.22±0.01 


0.21±0.01 


0.21±0.01 


A3112 


1.5 


2.9 


0.43±0.01 


0.49±0.01 


0.49±0.01 


0.53±0.01 


A3571 


0.0 


2.1 


1.95±0.02 


2.07±0.02 


1.93±0.02 


2.27±0.03 


A85 


1.5 


3.0 


1.10±0.02 


1.25±0.02 


1.18±0.02 


1.27±0.02 


Coma 


1.0 


5.0 


3.83±0.02 


4.29±0.02 


4.23±0.02 


4.34±0.03 


HydiaA 


1.5 


2.7 


0.52±0.01 


0.59±0.01 


0.54±0.01 


0.51±0.01 


MKW3S 


1.5 


2.5 


0.38±0.01 


0.37±0.10 


0.38±0.10 


0.42±0.01 



Notes. {a) The flux is given in units of 10~" erg s _I cnT 2 . The fluxes 
are scaled to correspond to the full annulus given by the inner and outer 
radius r in and r out . 



12.1. pn v.s. MOS 

We found that while the MOS2 hard band fluxes are in average 
only ~3% lower than those of MOS1, the difference is very sig- 
nificant (6.2<x, see Fig.[l6]and TablefTOt. There are bigger differ- 
ences between the pn and the MOS: MOS1 hard band fluxes are 
higher by ~7% (15.4<x) while MOS2 hard band fluxes are higher 



by 5~% (10. Oct) than those measured by the pn. The differences 
between the pn and MOS are similar in the soft band. This flux 
difference between the pn and MOS was reported earlier by e.g. 
Stuhlinger et al. (2008) and Mateos et al. (2009). While quali- 
tatively consistent with the analysis of the 2XMM catalogue in 
Mateos et al. (2009) (when considering similar off-axis angles 
(0-2 arcmin) and energy bands (0.5-12 keV)) our pn/MOS dif- 
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lines show the weighted mean of fp + the error of the mean. 



ferences are somewhat smaller. Our results are consistent with 
other IACHEC work based on a sample of blazars (Smith et al., 
2010, in prep.) 

12.2. EPIC v.s. ACIS 

We then compared the EPIC fluxes to those obtained using ACIS 
data and CALDB version 4.2. We found that the scatter of the 
flux between the pn and ACIS is larger (~6% ) around the mean 
value in both bands is much larger than the statistical uncertain- 
ties of the flux measurements (~1% ). This implies a systematic 
uncertainty component which varies from cluster-to-cluster. We 
examined this by restricting the spectral extraction regions to the 
common sky area where all the cameras have full coverage but 
this did not decrease the scatter significantly. We also experi- 
mented using two-temperature models with similar results. We 
further examined whether the differences in the shapes of the 
best-fit models between different instruments contribute to the 
flux scatter. We fitted the MOS and ACIS data with the best-fit 
pn model, allowing only the normalisation to be a free parame- 
ter. This did not reduce the scatter either. 

In the hard band, the ACIS fluxes are systematically and sig- 
nificantly higher than the pn fluxes, by 1 1+0.5% (24.7<x) on av- 
erage (see Fig. [PTl i. There is a trend of increasing relative flux 
difference with increasing temperature. ACIS flux values also 
exceed those of MOS1 and MOS2, but by a smaller amount: 
3±0.5% (6.1o-) and 6±0.5% (12.8cr), respectively. Thus, our 
analysis shows that the hard band fluxes using the pn, MOS or 
ACIS instruments, might differ by 5-10% due to similar uncer- 
tainty of the calibration of the normalisation of the effective areas 
in the hard band. 

In the soft band, the agreement between ACIS and the pn is 
much better than in the hard band: the fluxes differ only by ~2% 
on average (see Fig. [17}. This is interesting when noting that 
the soft band temperatures between EPIC and ACIS disagreed 
substantially (see Sect. [9]). 




pn hard band T [keV] pn hard band T [keV] pn hard band T [keV] 

Fig. 18. The upper panels show the best-fit temperatures (aster- 
isks) in the hard band obtained with XMM-Newton pn v.s. those 
obtained with Chandra ACIS using CALDB 3.4 (left panel), 
CALDB 4.1.1 (middle panel) and CALDB 4.2 (right panel). The 
lower panels show the relative ACIS-pn temperature difference, 
fj - AT/<T> (diamonds) and its lcr uncertainty. The dotted and 
dashed lines show the weighted mean of fp + the error of the 
mean. 



13. Chandra calibration changes 

In last few years (2007-2009) a considerable effort has been de- 
voted into re-analysing the ground-based calibration of Chandra 
. Here we examined how the temperatures and fluxes depend on 
the version of the Chandra calibration data base (CALDB), i.e. 
CALDB 3.4 (Dec 2006), CALDB 4.1.1 (Jan 2009) and CALDB 
4.2 (Dec 2009), (see e.g. David et al., 2007; David et al., 2009). 

13.0.1. CALDB 3.4 

After the Chandra launch, gratings observations of AGN showed 
that the mirror reflectivity near the Ir-M edge was 15% higher 
than predicted by the HRMA effective area model. In CALDB 
3.4, a new HRMA effective area model was released which in- 
cluded a 22 A layer of hydrocarbon contaminant on the mir- 
rors and reproduced the observed reflectivity across the Ir-M 
edge. In our earlier work (Nevalainen et al. 2007) we found that 
this model yielded higher hard band temperatures for clusters 
hotter than ~4 keV when compared to values derived with the 
XMM-Newton EPIC detectors using the calibration as of Aug 
2007. Similar results have been obtained with an older version 
of Chandra calibration (Vikhlinin et al., 2005; Kotov et al., 2005; 
Snowden et al., 2008). 

In the current work we found that the Chandra CALDB 3.4 
hard band temperatures exceeded those obtained with the pn us- 
ing the latest EPIC calibration information in Dec 2009 by an 
average of ~10% (see Figs.[18]and[T9]and Table lA.ll . The max- 
imum temperature difference between the EPIC Dec 2009 cali- 
bration and ACIS CALDB 3.4 is -30%. 

13.0.2. CALDB 4.1.1 

Following our cluster temperature Chandra /XMM-Newton 
cross-calibration work (Nevalainen et al. 2007), the ground- 
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Fig. 19. The relative differences of the soft and hard band temperatures (two upper rows) and fluxes (two lower rows) between the 
pn and ACIS at different stages of ACIS calibration (different columns). The values are shown as a function of the hard band pn 
temperatures. Note that the scaling of the y-axis is different in different rows. The average values are marked with a dotted line. 



based HRMA calibration data taken at the X-ray Calibration 
Facility (XRCF) at MSFC was re-analysed. It was found that 
the molecular contamination was already present on the mirrors 
at XRCF and that the correction applied to the CALDB3.4 ver- 
sion of the HRMA effective area over-corrected the effect of the 
contamination. This problem was corrected in CALDB 4.1.1. 
During XRCF testing, shutters were used to calibrate each of 
the four mirror shells independently. CALDB 4.1.1 contains a 
model with the XRCF measured depths for the contaminant on 
each shell. 

As a consequence, the change in the reflectivity above 2 keV 
produced lower hard band temperatures. We found in the cur- 
rent work that compared to XMM-Newton pn values, the ACIS 
temperatures obtained with the CALDB 4.1.1 agreed within a 
few per cent (see Figs. [T8land[T9l. Also, the soft and hard band 
fluxes increased by ~9% and ~4%, respectively. 

13.0.3. CALDB 4.2.0 

ACIS observations of astronomical sources as well as its own 
external calibration source (ECS) have shown that molecular 
contamination has been building up on the optical/UV blocking 
filters since launch. A new version (N0005) of the ACIS con- 
tamination model was released in CALDB 4.2.0. The previous 
version of the ACIS contamination model provided a good fit 



to the ECS data up until about 2006. For more recent observa- 
tions, the previous version underestimated the optical depth of 
the contaminant on the ACIS filters. The new version released 
in CALDB 4.2 produces a much better agreement with the ECS 
data for recent observations. 

We found that this revision mostly affected the cluster soft 
band fluxes, and produced consistency with the pn, while the 
soft band temperatures did not change appreciably (see Fig.[T9b 

14. Conclusions 

We performed an X-ray spectral analysis of clusters of galax- 
ies using data obtained with the XMM-Newton EPIC instruments 
pn, MOS1 and MOS2 and Chandra ACIS-S and ACIS-I. We ad- 
ditionally used the published results obtained using BeppoSAX 
MECS (deGrandi & Molendi, 2002). We compared the derived 
temperatures and fluxes for each cluster based on the current in- 
strument calibrations as of December 2009. The results are sum- 
marised in Fig. [20] and Table [TT] 

We found that there are no systematic differences in the tem- 
peratures obtained by fitting the 2-7 keV energy band of ACIS 
and EPIC (2-10 keV for MECS). These values are also consis- 
tent with those obtained by the EPIC temperature measurements 
based on the Fe XXV/XXVI line ratio. This shows that the en- 
ergy dependence of the effective area in this band is accurately 
modelled in the studied instruments. Thus, the IACHEC cluster 



17 



J. Nevalainen et al.: IACHEC cluster calibration 



A 
h— 
V 



< 

Li_ 

V 



0.2 r 



0.1 r 



< -0.1 r 



-0.2 r 



T 1 1 1 — 

SOFT BAND <> 



❖ 



M1/PN M2/PN M2/M1 ACIS/PN 



i 1 1 1 r 

HARD BAND 



o.o r • ~i-~~r-- 



<5> 



M1/PN M2/PN M2/M1 ACIS/PN MECS/PN J 



AF/<F> 

AT/<T> 



J I I I I I I I I L 



Fig. 20. The average relative difference (diamonds) + the error of the mean of the fluxes (solid line) and temperatures (dotted line) 
for different instrument pairs in the soft band (left side of the plot) and in the hard band (right side of the plot). 

Table 11. The relative difference and its significance of temperatures and fluxes for different combinations of instruments and 
methods. 
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Notes. The values correspond to the best-fit single temperature models in the <a) 2-7 keV band, 0.5-2 keV band, (c) 0.5-7 keV band, (rf) ~6. 0-7.0 
keV band and <e) 2-6 keV band The ACIS values are obtained using CALDB 4.2.0. <e) ji gives the weighted mean of the relative difference between 
the two measurements in percentages. tf) sig. gives the significance of the relative difference in terms of cr. 



sample in the studied radial range (~ 0.1-0.3 rgoo) constitutes a 
set of standard candles for the calibration of the energy depen- 
dence of the hard band effective area. On the other hand, the 
disagreements at 6-25<r level on the hard band fluxes showed 
that there are systematic calibration uncertainties in the normal- 
isations of the effective areas by 5-10% in the 2-7 keV band. 

The temperatures obtained by fitting the soft band (0.5-2.0 
keV) data of the pn and ACIS differ systematically and signif- 



icantly, by ~18% (i.e. at 8.6<r level) on average. This indicates 
remaining uncertainties with the calibration of the energy de- 
pendence of the effective area in the 0.5-2.0 keV band in one 
or all instruments. Comparison of the residuals showed that the 
relative pn/ACIS cross-calibration bias is approximately a linear 
function of energy in the 0.5-2.0 keV band, amounting to a vari- 
ation of 10% in this band. The uncertainties of the calibration 
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and the modelling of the cluster emission in the soft band cannot 
be simultaneously resolved by using clusters of galaxies alone. 

Due to the higher effective area and the higher number of 
intrinsic cluster photons in the soft band, the statistical weight 
of the soft band data is much higher than that of the hard band. 
Thus, the calibration uncertainties in the soft band will affect 
the scientific analysis of clusters of galaxies, when using the full 
useful energy band (~0.5-7.0 keV). Considering the variation 
between the different instruments, the 0.5-7.0 keV band tem- 
perature measurement of clusters of galaxies with EPIC/XMM- 
Newton or AClS/Chandra is uncertain by 10-15% on average. 
These uncertainties will also affect the analysis of the wide band 
continuum spectra of other types of objects using ACIS or EPIC. 

We evaluated the systematic effects on the Fe XXV/XXVI 
line ratio temperature measurement due to uncertainties of the 
implemented EPIC calibration of the energy scale and energy 
resolution and redistribution and the details of the line ratio 
modeling. The effect on the measured temperature is ~4%. 
The temperatures measured using the continuum shape and the 
Fe XXV/XXVI line ratio agree very accurately. This indicates 
that the deviations from the ionisation equilibrium state and 
Maxwellian electron velocity distribution are negligible in the 
studied regions of this cluster sample. Since the Fe XXV/XXVI 
line ratio measurement is not sensitive to smooth changes in the 
calibration of the effective area, it could be a powerful hard band 
calibration tool in future X-ray missions. 

After the submission of this paper an update on the pn redis- 
tribution become public (Haberl et al., 2010). This may improve 
the Fe XXV/XXVI emission line modeling and we will address 
this issue in a follow-up paper. The effect on the fluxes and con- 
tinuum temperatures is likely small, because the line emission 
(besides Fe XXV and XXVI) is weak in the cluster sample. A 
study of a XMM-Newton sample of sources with a wide range 
of spectral shapes (Stuhlinger et al. 2010, in prep.) shows that 
the change of the relative pn/MOS flux above 0.5 keV between 
the new calibration and that used in the present paper is 3% at 
the most. Thus, the cluster flux discrepancy between the pn and 
MOS will not be resolved with the new pn calibration in July 
2010. 
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Appendix A: Spectral fits 

We show here the data and the best-fit single-temperature 
MEKAL models in the hard and soft bands for pn (Figs. IA.1I 
andlAT2li and for ACIS (Figs.lA31andlA4li used for the pn/ACIS 
temperature and flux comparison. The spectral parameters of the 
pn and ACIS fits in all bands are shown in Table I A. fl 
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normalised counts s 1 keV 




Fig. A.l. The hard band pn spectra (crosses) and the best-fit single-temperature fits (solid lines) for the cluster sample. The normal- 
isation of the spectra is adjusted for plot clarity and does not correctly reflect the relative brightness of the clusters. 
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Fig. A.2. The soft band pn spectra (crosses) and the best-fit single-temperature fits (solid lines) for the cluster sample. The normali- 
sation of the spectra is adjusted for plot clarity and does not correctly reflect the relative brightness of the clusters. 
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Fig.A.3. The hard band ACIS spectral data (crosses) and the best-fit single-temperature fits (lines) for the cluster sample. The 
normalisation of the spectra is adjusted for plot clarity and does not reflect correctly the relative brightness of the clusters. The 
spectra for A85 and Coma were obtained from several CCD chips and were fitted simultaneously. 
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Fig.A.4. The soft band ACIS spectra (crosses) and the best-fit single-temperature fits (solid lines) for the cluster sample. The 
normalisation of the spectra is adjusted for plot clarity and does not correctly reflect the relative brightness of the clusters. The 
spectra for A85 and Coma were obtained from several CCD chips and were fitted simultaneously. 



24 



J. Nevalainen et al.: IACHEC cluster calibration 
Table A.l. Best-fit temperatures and metal abundances for pn and ACIS in different bands 
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101.6/99 



SOFT BAND (0.5-2.0 keV) 



A 1795 

A2029 

A2052 

A2199 

A3112 

A3571 

A85 

Coma 

HydraA 

MKW3S 



1.5 
1.5 
1.7 
2.0 
1.5 
0.0 
1.5 
1.0 
1.5 
1.5 



2.7 
2.5 
2.5 
2.9 
2.9 
2.1 
3.0 
5.0 
2.7 
2.5 



4.7[4.5-4.8] 
6.1 [5.7-6.6] 
2.7[2.6-2.8] 
3.8[3.7^1.0] 
3.9[3.7-4.1] 
5.7[5.3-6.0] 
5.2[4.9-5.6] 
6.0[5.8-6.1] 
2.9[2.7-3.0] 
3.2[3. 1-3.3] 



0.33[0.29-0.36] 
0.30[0.19-0.41] 
0.48 [0.45-0.52] 
0.40[0.35-0.45] 
0.35[0.29-0.41] 
0.39[0.3 1-0.47] 
0.50[0.41-0.60] 
0.19[0.15-0.22] 
0.28 [0.24-0.33] 
0.36[0.33-0.40] 



232.6/303 
203.1/231 
226.7/250 
251.8/255 
196.5/218 
258.3/276 
208.3/220 
325.7/303 
106.7/131 
220.6/277 



0[5.6-6.5] 
0[7.7-8.3] 
2[3.2-3.3] 
0[4.7-5.2] 
1 [4.6-5.6] 
7[7.0-8.5] 
9[4.7-5.1] 
9[8.5-9.4] 
3[3.3-3.4] 
3[3. 1-3.4] 



0.34[0. 
0.3 1[0. 
0.61 [0. 
0.70[0 
0.33[0 
0.69[0. 
0.23 [0 
0.63[0. 
0.37[0. 
0.32[0. 



24-0.45] 
23-0.40] 
58-0.64] 
63-0.80] 
21-0.46] 
51-0.91] 
17-0.28] 
49-0.82] 
35-0.40] 
28-0.36] 



115.6/100 
115.7/100 
165.8/100 
104.9/100 

106.1/82 
110.2/100 
200.8/180 
368.8/315 

96.7/100 
116.3/100 



WIDE BAND (0.5-7.0 keV) 



A 1795 

A2029 

A2052 

A2199 

A3112 

A3571 

A85 

Coma 

HydraA 

MKW3S 



1.5 
1.5 
1.7 
2.0 
1.5 
0.0 
1.5 
1.0 
1.5 
1.5 



2.7 
2.5 
2.5 
2.9 
2.9 
2.1 
3.0 
5.0 
2.7 
2.5 



5.4[5.4-5.5] 
7.0[6.8-7.2] 
3.0[2.9-3.0] 
4.2[4. 1-4.3] 
4.4[4.3-4.5] 
6.6[6.4-6.7] 
5. 5 [5 .4-5 .7] 
7.1[7.0-7.2] 
3. 3 [3. 2-3 .4] 
3.6[3.5-3.6] 



0.44[0.42-0.47] 
0.47[0.42-0.53] 
0.56[0.53-0.60] 
0.49[0.45-0.53] 
0.41 [0.36-0.45] 
0.53 [0.49-0.57] 
0.53[0.47-0.58] 
0.28 [0.27-0.30] 
0.33[0.28-0.39] 
0.42[0.41-0.45] 



509.5/585 
294.4/345 
308.2/343 
368.4/368 
270.7/308 
395.7/444 
272.0/317 
969.4/862 
152.9/181 
335.6/413 



3 [6.2-6.5] 
9[8.8-9.1] 
3[3.3-3.4] 
7[4.6-4.8] 
7[5.5-5.9] 
6[8.4-8.8] 
8[5.7-5.9] 
0[8.8-9.1] 
8[3.7-3.8] 
8[3.7-3.8] 



0.42[0. 
0.52[0. 
0.60[0. 
0.57[0. 
0.65[0. 
0.72[0. 
0.46[0. 
0.40[0. 
0.41 [0 
0.44[0. 



37-0.49] 
48-0.55] 
57-0.62] 
52-0.62] 
55-0.76] 
64-0.81] 
43-0.51] 
36-0.44] 
39-0.43] 
41-0.48] 



208.4/188 
398.1/360 
387.5/295 
200.4/181 
180.7/120 
176.7/191 
415.3/347 
830.2/753 
312.9/289 
241.1/201 



Notes. Temperatures and their uncertainties at lcr level were obtained with a single-temperature MEKAL model using data from the pn and ACIS 
instruments in different bands 
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Appendix B: Systematic uncertainties affecting the 
Fe XXV/XXVI line ratio measurement using 

XMM-Newton EPIC instruments 

B.1. Effective area 

B.1 .1 . Break in the effective area 

The virtue of the Fe XXV/XXVI line ratio method is that the 
energies of the two lines only differ by ~300 eV. If there was 
a sudden change in the real effective area, co-incidentally at an 
energy between those of the FeXXV and FeXXVI line energies, 
that was not correctly implemented into the calibration, then the 
measured line ratio would be inaccurate. However, the redshifts 
in our sample vary from 0.023 1 to 0.0773 which produce a ~300 
eV variation in the Fe line energies in different clusters. Thus, a 
possible break in the effective area at a fixed energy cannot pro- 
duce a systematic effect in the sample, but can affect some of 
the clusters. To estimate the effect of this unlikely situation, we 
simulated a pn spectrum of a cluster with T=9 keV and metal 
abundance of 0.3 Solar at a redshift z=0. When fitting the sim- 
ulated data in a 6.45-7.25 keV band, we introduced a 10% drop 
in the effective area in the associated auxiliary response file at 
E>6.85 keV. In order to produce the data above 6.85 keV, the 
model Fe XXVI flux increases, i.e. the temperature increased to 
10.4 keV, which is 15% higher than the input value. 

B.1.2. Smooth calibration bias 

We note that we found an excellent agreement between the hard 
band temperatures obtained with different instruments (see Sect. 
[7]). Because the different instruments are calibrated based on 
a combination of ground measurements and in-flight measure- 
ments of different types of celestial sources, it would be a highly 
unlikely co-incidence that all the instruments have similar un- 
certainties in the calibration as to yield similarly biased temper- 
atures. This hypothetical bias should be a smooth function of 
energy in order for the biased model to agree well with the data, 
as observed. We examined the sensitivity of the Fe XXV/XXVI 
line ratio temperature measurement to a smooth hard band cali- 
bration bias using simulated data. We created spectra with kT in 
the range 4-10 keV. We kept the metal abundance at 0.3 Solar 
and the redshift at 0. We used 10 7 counts in the 6.45-7.25 keV 
energy band in order to eliminate the statistical effects. We used 
the auxiliary response file of A1795 for the simulations. 

When fitting the simulated data, we approximated the effect 
of a smooth calibration bias by modifying the auxiliary file we 
used for the simulations. In detail, we multiplied the effective 
area column with a linear function which has a value of 1 .0 at 1 
keV, and varies between 0.8 and 1.2 at 10 keV. Thus the effec- 
tive area is unchanged at 1 keV, but is underestimated or overes- 
timated by a maximum of 20% at 10 keV. 

We found that the resulting relative bias in the temperature 
measurement is higher for the higher input temperatures, due 
to their higher statistical weight in the more biased high energy 
band (see Fig. |BT]). However, the effect in the Fe XXV/XXVI 
temperature measurement was very small: even if the effective 
area was off by 15% at photon energies of 10 keV, the best-fit Fe 
XXV/XXVI temperature would be biased by less than 1 % (see 
Fig. IB.U . The reason for this is that even though the effective 
area is over- or underpredicted by up to ~ 15% in the 6-7 keV 
band in average, the relative change of the effective area between 
the 6.45 and 7.25 keV is only ~1%. Thus the only effect is a 
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Fig. B.l. The temperature bias in the simulated spectra as a func- 
tion of the fraction of the linear bias in the effective area at 10 
keV for the 2.0 - 7.0 keV band fits (dotted lines) and for the ~6- 
7 keV band fit (dashed lines). The different dotted lines denote 
different input temperatures from 4 to 10 keV (higher temper- 
atures yield higher bias). For the ~6-7 keV band fit only the 
curves for kT=10 keV are shown for clarity. The lines in the up- 
per half (lower half) correspond to positive (negative) bias in the 
effective area. 

biased MEKAL normalisation up to ± ~10%, while the shape of 
the spectrum does not change. 

The situation is the opposite in the 2.0-7.0 keV band. The 
varying degree of calibration bias in a 10-20% range does not 
much affect the derived emission measures (which only vary by 
~1%). However, it does affect the measured continuum shape, 
yielding a 10-20% bias in the temperature (see Fig. IB. II ). 

B.1.3. ACIS calibration changes 

We further tested the effects of smooth effective area calibration 
changes on the Fe XXV/XXVI line ratio measurement by a com- 
parison of ACIS temperatures derived using calibration versions 
CALDB 3.4 and CALDB 4.2.0. Note that the substantial, smooth 
changes in the ACIS effective area calibration between these two 
versions resulted in significant changes in the hard and soft band 
temperatures (see Sect.[7]and[9]). We relaxed the requirement for 
the minimum number of counts to 400 since the bias due to the 
limited number of counts (see Sect. I8.11 l does not depend on the 
calibration version and thus will not introduce a systematic dif- 
ference. 

We fitted the data in the ~6-7 keV energy band allowing 
the emission measure to be a free parameter. We found that the 
Fe XXV/XXVI best-fit temperatures obtained with CALDB 3.4 
and CALDB 4.2 are nearly identical (see Fig. IB.2t . This result 
supports the above suggestion that the Fe XXV/XXVI line ratio 
temperature is rather insensitive to smooth changes in the cali- 
bration of the effective area. 

B.2. Accuracy of the energy resolution and redistribution 
calibration 

We examined how robust our Fe XXV/XXVI temperature mea- 
surements are when considering possible uncertainties in mod- 
elling of the time-dependent energy resolution and redistribu- 
tion of the pn and MOS (Guainazzi, 2010). Following the work 
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Fig.B.2. The best-fit temperatures obtained from ACIS data in the ~6-7 keV energy band using calibration versions CALDB 3.4 
v.s CALDB4.2.0. 



of Molendi et al. (2009), we added a multiplicative "gsmooth" 
component to the model, which convolves the input model with 
a Gaussian kernel. We first fitted each cluster individually. This 
yielded rather large statistical uncertainties for the width of the 
Gaussian kernel (= cr) within which the clusters are rather con- 
sistent. On average, the value for cr is 26 eV, and the rms scatter 
is 20 eV. Thus, within the scatter, these results are roughly con- 
sistent with the Gaussian width of 4e V reported by Molendi et 
al. (2009) for Perseus. 

However, the cluster observations of our sample were per- 
formed in a short period, early in the XMM-Newton mission (dur- 
ing years 2000-2002) when the energy resolution was stable (e.g. 
Guainazzi 2010). We thus experimented by simultaneously fit- 
ting the spectra of all clusters including the gsmooth-component, 
while forcing cr to be equal in all data sets. For x 1 - compar- 
ison, we also repeated the fit without the gsmooth-component. 
We found that the addition of the gsmooth-component improved 
the fit significantly: x 1 decreased from 297.0 to 285.9 while the 
number of free parameters increased from 18 to 19 while us- 
ing 345 spectral bins. The fit obtains a value of 35+7 eV for 
the Gaussian width, i.e. consistent with that obtained when fit- 
ting each cluster individually above. This value is consistent 
with measurements of the on-board pn calibration source (see, 
e.g., Guainazzi et al. 2010). Thus, our results indicate a system- 
atic underestimation of the pn resolution in the response matrix 
used in this paper. A similar exercise to our MOS data yielded 
a Gaussian width consistent with zero, i.e. this effect is not ap- 
parent in MOS (consistent with Guainazzi et al., 2010). Thus 
this is not a physical broadening effect of the lines, but rather an 
instrumental problem of the pn. 

The pn temperatures decreased systematically with the in- 
clusion of the gsmooth-component, by 2% on average. This is 



slightly smaller than the value obtained when fitting each cluster 
individually. Thus we conclude that the effect of the uncertain- 
ties in the calibration of the energy resolution and redistribution 
at ~ 6-7 keV render the measured Fe XXV/XXVI pn tempera- 
tures too high by ~2%. 



B.3. Calibration of the energy scale 

We then examined how possible problems in the calibration 
of the energy scale in the ~6-7 keV band impact the Fe 
XXV/XXVI temperature measurements. There are two principal 
factors affecting the energy reconstruction of an event, namely 
Charge Transfer Inefficiency (CTI) and gain. CTI describes the 
transfer of charge as it is transported through the CCD to the 
output amplifiers during read-out. Gain is the conversion of the 
charge signal deposited by a detected photon from charge into 
energy (Guainazzi 2010). Note that possible uncertainties in the 
redshift measurement using the spectra of the galaxies of a given 
cluster (NASA Extragalactic Database) may contribute to the un- 
certainty of the energy scale because we keep the redshift fixed 
in the fits. 

In order to examine the above effects, we used the "gain 
fit" option in XSPEC 11.3.2.ag package to modify the defini- 
tion of the channel energies in the response matrix with an off- 
set that does not depend on the energy. We fitted each cluster 
individually and found that in average the offset parameter ob- 
tained values of 5 eV and 16 for the pn and MOS, respectively. 
Considering the rms scatter, these values are consistent with the 
nominal accuracies of energy reconstruction of 10 eV and 5 
eV for the pn and MOS, respectively, in the full energy range 
(Guainazzi 2010). 
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The modification of the gain offset renders the temperatures 
systematically lower, by 1% (3%) for the pn (MOS) on average. 
Thus, we conclude that the uncertainties in the calibration of the 
energy scale of the EPIC instruments combined with the possi- 
ble uncertainties of the cluster redshift measurements render the 
measured Fe XXV/XXVI EPIC temperatures too high by ~1- 
3%. 

6.4. Details of the emission model 

In addition to the MEKAL model we adopted for the emis- 
sion modelling, there also exists a commonly used emission 
model APEC (Smith et al., 2001). Differences on the details 
of the modelling of electron transitions may affect the inter- 
preted temperature for a given line ratio data. We examined 
this by fitting the ~6-7 keV band pn data using either MEKAL 
or APEC model. We found that the temperatures derived with 
APEC model are systematically lower than those obtained with 
the MEKAL model, by 2% on average. 

We also experimented with the choise of the element abun- 
dance tables by replacing our adopted one(Grevesse & Sauval, 
1998) with that of Anders & Grevesse (1989) or Lodders (2003). 
The differences in the derived ~6-7 keV band temperatures are 
smaller than 0.1%. 
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